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ABSTRACT 

A  DYNamic  Tuibinc  Engine  Compressor  Code  (DYNTECC)  has 
been  modified  to  model  speed  transients  fiom  0-100%  of  ocmipressor 
design  speed.  Hie  impetus  for  this  enhancement  was  to  investigate 
stage  matching  and  stalling  behavior  during  a  start  sequence  as 
axnpared  to  rotating  stall  events  above  ground  idle.  The  model  can 
simulate  speed  and  throttle  excursions  simultaneously  as  well  as  time 
varying  bleed  flow  schedules.  Results  of  a  start  simulation  are 
presented  and  compared  to  experimental  data  obtained  from  an  axi- 
centrifugal  turboshaft  engine  and  companion  compressor  rig.  Stage 
by  stage  comparisons  reveal  the  fitwit  stages  to  be  operating  in  or  near 
rotating  stall  through  most  of  the  start  sequence.  The  model  matches 
the  starting  operating  line  quite  well  in  the  forward  stages  with 
deviations  qipearing  in  the  rearward  stages  near  the  start  bleed. 
Overall,  the  performance  of  the  model  is  very  promising  and  adds 
significantly  to  the  dynamic  simulation  capabilities  of  DYNTECC. 


NOMENCLATURE 

A  Flow  Area 
e  Internal  Energy 
F  Blade  Force 
H  Total  Enthalpy 
M  Mach  Number 
P  Pressure 
Q  Heat  Addition 
R  Ideal  Gas  Constant 
S  Shaft  Work 
T  Temperature 
t  Time 

U  Axial  Flow  Velocity 
W  Mass  Flow  Rate 
X  Axial  Coordinate 
p  Density 
y  Specific  Heat  Ratio 


Subscripts 

B  Bleed  Flow 

8  Static  Property 

t  Total  (Stagnation)  Property 

X  Axial  Direction 


INTRODUCTION 

The  development  of  reliable  high  performance  gas  turbine 
engines  continues  to  be  one  of  the  most  challenging  engineering 
endeavors  the  20^^  century.  In  today's  competitive  market,  engine 
designers  rely  increasingly  on  the  ability  to  numerically  simulate 
engine  performance  throughout  the  entire  region  of  operation.  To  this 
end,  considerable  cff(xt  has  been  devoted  to  modeling  the 
compression  system  ol  gas  turbines  with  particular  attention  to 
multistage  axial  compressors  for  aircraft  applications.  The  work 
presented  here  stems  fiom  a  consortium  of  government,  industrial, 
and  academic  members  known  as  the  Joint  Dynamic  Airbrcalhing 
Propulsion  Simulations  (JDAPS)  partnership  whose  mission  is  to 
advance  the  state  of  the  art  in  numerical  modeling  of  gas  turbine 
engine  components  (Davis  et  al.,  1995). 


BACKGROUND 

Numerous  publications  in  the  past  two  decades  have  been 
devoted  to  understanding  the  phenomenon  of  surge  and  rotating  stall 
in  axial  flow  multistage  compressors.  Numerical  models  capable  of 
predicting  post  stall  behavior  generally  gained  acceptance  with  the 
lumped  volume  ^jproach  of  Greitzer  (1976)  $nd  have  evolved  into 
finite  difference  control  volume  methods  which  can  isolate 
aerodynamic  behavior  of  an  individual  stage.  One  widely  accepted 
and  validated  model  known  as  DYNTECC  was  developed  by  Davis 
and  O’Brien  (1986,1991)  and  has  become  the  cornerstone  of  the 
JDAPS  modeling  effort 

An  application  of  this  model  by  O’Brien  and  Boyer  (1989)  was 
able  to  identify  the  critical  stages  in  a  10  stage  high  performance 
compressor  which  exhibited  a  rotating  stall  pit^lcm  at  mid  speeds. 
One  conclusion  of  their  results  was  that  stall  recovery  problems  at 
mid-range  speeds  can  be  identified  as  an  extension  of  the  starting 
problems  typically  encountered  fiom  zero  speed  to  ground  idle.  This 
"extended  starting’*  theory  postulated  that  methods  which  provide  for 
aerodynamic  starting  at  low  speeds  will  produce  recovery  at  higher 


spe  eds  if  pn^ly  applied.  The  woik  presented  here  was  originally 
undertaken  to  investigate  this  extended  starting  theory  from  Ae 
perspective  of  a  <fynamic  model.  The  T55-L-712  turboshaft  cngiiw 
was  chosen  as  the  validation  vehicle  for  this  model  because  a 
timely  start  testing  program  begun  in  1993  which  included  interstage 
data  frOTi  20%’100%  of  design  speed  (Owen  and  Davis,  1994).  The 
availabiliQr  oi  complete  interstage  data  in  the  starting  region  was 
crucial  to  die  model  development 

The  investigation  of  starting  problems  and  the  development  of 
full  speed  range  models  has  been  given  some  attention  in  current 
literature,  CSiappcU  and  l^Laughlin  (1993)  developed  m  ^ipro^ 
for  modeling  continuous  turbine  engine  operation  whidi  was  an 
extension  of  a  component  matching  technique  already  in  use. 
Agrawal  and  Yunis  (1982)  presented  a  similar  approadi  based  on 
generalized  component  mqw  in  the  starting  region.  These  models 
generally  depend  on  table  lookiq)  or  “mq)  reading”  schemes  and  do 
not  have  the  inherent  ibiliiy  to  lode  at  post  stall  and  interstage 

behavior.  Some  Navicr-Stokes  calculations  have  been  undertaken  for 

starting  flows  in  a  compressor  cascade  (Outa  ct  al.,  1993)  but  the 
computing  time  and  memory  requirements  presently  make  this 
method  impractical  for  modeling  a  multistage  compressor  with  post 
stall  behavior.  The  utility  of  the  DYNTECC  starting  model  lies  in  its 
ability  to  look  at  c<wnprcssor  behavior  for  any  foreseeable  operating 
condition  based  on  fundamental  flow  j^ysics  with  reasonable 
computing  lime. 


STARTING  MODEL 

DYNTECC  is  a  one-dimensional  stage-by-stage  compression 
system  model  which  is  able  to  analyze  generic  compression  systems 
as  illustrated  in  Rg.  1.  The  inlet  ducting,  compressor,  and  cewnbustor 
vdumes  are  modeled  as  a  series  erf  elemental  control  volumes  with 
mass,  woA,  and  heat  fluxes  across  the  control  volume  surfaces.  The 
equations  governing  the  flow  are  those  for  mass,  momentum,  and 
energy  transfer  for  a  non-viscous  fluid  with  turbomachincry  source 
terms,  commonly  referred  to  as  the  Euler  equations  (Eq.l). 

The  source  terms  for  cempressor  stages  are  given  in  the  form  of 
pressure  and  temperature  coefficients  vs.  mass  flow  coefficient  for 
each  stage.  The  stage  characteristics  are  defined  for  all  performance 
possibilities,  which  arc  generally  divided  into  three  distinct  regions 
as  shown  in  Fig.  2.  The  pre-stall  region  is  the  normal  stable 
operating  regime  where  the  pressure  rise  characteristic  has  a  negative 
slope.  The  rotating  stall  region  is  modeled  as  a  continuous 
characteristic  along  a  throttle  line  with  positive  slope  down  to  the 
zero  flow  condition,  and  the  reversed  flow  region  represents 
performance  associated  with  full  annulus  reversed  flow.  The  results 
presented  here  arc  concerned  primarily  with  normal  operation  up  to 
the  inception  of  rotating  stall.  Ibis  is  the  region  in  which 
experimental  stage  characteristics  were  available  from  the  T-55  test 
program.  A  full  complcm<mt  of  spanwisc  pressure  and  temperature 
pro^  for  each  stage  was  used  to  measure  stage  performance  from 
20%A00%  erf  design  cometed  speed,  with  ground  idle  occurring  at 
approximately  61%.  It  should  be  noted  that  the  flow  blockage  due  to 
boundary  layer  effects  will  be  inherent  in  the  stage  characteristics, 
thus  blo^age  effects  are  not  handled  explicitly  inside  the  framework 
of  DYNTECC. 


a)  Compressor  and  ducting  system 
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Figure  1.  DYNTECC  Control  Volume  Technique. 
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Rgure  2.  Typical  Stage  Characteristics 
BOUNDARY  CONDITIONS 

Some  changes  to  the  model  have  been  implemented  in  order  to 
simulate  a  start  sequence  from  0%  -  100%  design  corrected  speed. 
The  most  critical  these  was  the  development  of  a  xmiform 
pressurc/arca  boimdary  condition  which  can  be  applied  to  both 
choked  and  unchoked  exit  flow.  This  allowed  the  model  to  nm  at 
variable  speeds  over  the  entire  range  of  operability.  The  same 
principle  was  applied  to  develop  a  start  bleed  model  based  on  bleed 
flow  area  and  bleed  dump  pressure.  This  allowed  simulation  of 
actual  start  bleed  schedules  and  significantly  enhanced  the  ability  to 
model  interstage  behavior  during  a  start 


The  inflow  boundary  condition  during  normal  forward  flow  is 
the  specification  of  total  pressure  and  temperature.  The  exit  boundary 
condition  has  traditionally  been  the  ^xjcification  of  exit  mass  flow 
parameter  or  static  pressure.  The  ^jccification  of  exit  mass  flow 
parameter  makes  the  implicit  assumption  of  a  choked  downstream 
throttle  area,  while  q)ecification  of  static  pressure  assumes  unchoked 
conditions.  This  set  of  boundary  conditions  works  well  at  constant 
speed  as  DYKTBOC  was  originally  developed,  but  proved  to  be 
cumbersome  when  varying  speeds  from  0%-100%  of  design  speed, 
since  it  is  unknown  a  priori  at  what  point  the  exit  fiow  will  become 
choked  as  speed  increases. 

To  model  starting  behavior  over  a  wide  range  of  speeds,  a 
uniform  set  of  boundary  conditions  was  developed  based  on  the  exit 
flow  area  and  exit  static  pressure.  This  provides  a  logical  model 
since  the  actual  flowfield  is  also  governed  by  these  two  physical 
properties.  The  pressure  /  area  boundary  condition  provides  all  the 
information  for  closure  of  the  governing  equations  and  ronatns  valid 
over  the  entire  operating  region.  The  exit  mass  flow  parameter  is 
still  employed  within  the  framewOTk  of  DYNTECXZ,  after  being 
derived  from  the  pressure  area  conditions  by  Eqs.  2  &  3.  The  exit 
Mach  number  is  determined  from  the  given  static  pressure  snd 
calculated  total  pressure  by  Eq.  2.  This  is  then  used  to  calculate  an 
exit  mass  flow  parameter,  which  incorporates  die  given  exit  flow  area 
in  Eq.  3.  This  technique  allows  for  speed  and  throttle  transients  over 
the  entire  compressor  performance  map  with  a  smooth  transition 
through  the  starting  region. 


START  BLEED  MODEL 

The  same  procedure  used  for  determining  the  exit  mass  flow  was 
also  used  to  develop  a  dynamic  bleed  flow  model.  Given  the  bleed 
flow  area  and  dump  static  pressure,  the  mass  flow  through  the  Weed 
is  determined  for  each  time  step.  For  the  T55  model,  the  ratio  of 
static  pressure  at  the  sixth  stage  start  bleed  to  the  bleed  dump  static 
pressure  determines  the  bleed  mass  flow  rate  for  a  given  bleed  area. 
This  allows  a  bleed  flow  schedule  to  be  modeled  based  on  the  actual 
bleed  flow  area  as  indicated  by  the  position  of  the  bleed  band  during 
a  start  sequence.  Bleed  flow  blockage  and  pressure  bsscs  through 
the  bleed  band  and  piping  arc  not  modeled  explicitly,  but  arc 
included  implicitly  in  the  selection  of  bleed  flow  area  and  bleed 
dump  pressure.  The  combination  of  exit  boundary  conditions  md 
bleed  flow  boundary  conditions  dclcnninc  the  internal  stage  matching 
in  the  model. 


ENGINE  DATA  ACQUISmON 

bi  1990,  the  U.S.  Anny  Vehicle  I^opuliion  Directonte  (VTO) 
initiated  its  non-ieoovcrable  stall  program.  One  goal  of  this  broad 
program  kx^dng  at  cfynamic  engine  events  is  the  stuefy  of  the 
AlliedSignal  T55*L*712  turboshaft  engine  and  to  report  to  the  U.S. 
Anny  Aviadon  and  Itoop  Command  (ATCOM)  on  this  engine’s  start 
sequence.  To  accomplish  this,  the  VPD’s  program  included 
extensive  rig  testing  which  de&ied  the  mdividual  stage 
characteristics  and  engine  testing  to  characterize  the  engine  start 
sequence. 

Engine  testing  began  in  November  of  1994  and  start  testing 
ended  on  March  27, 1995.  During  that  time,  127  engine  staru  were 
aoooaq>liahed  with  high  response  data  acquired  during  74  of  those 
starts.  Ihe  engine  was  started  at  altitudes  up  to  qjproximately  45 
bn,  irith  varying  fuel  schedule,  starter  torque  and  disengagement 
pennt,  turbine  inlet  temperature  (warm  restarts),  and  bleed  flow. 
Owen  (1995)  presents  in  detail  the  1^5-L-712  engine, 
instrumentation  locations,  daU  acquisition  procedures,  and  provides 
some  preliminaxy  results.  Data  from  this  program  is  available,  under 
appropriate  proprietary  constraints,  to  the  United  States 
turbonmehinery  community. 


transducers  located  on  the  shroud  in  front  of  the  first  rotor,  and  on  the 
hub  at  the  exit  of  die  first  two  stators. 

Study  state  pressure  is  sensed  using  an  electrically  scanned 
pressure  (ESP)  system  which  scans  and  iqxlates  the  readings  every 
second.  Start  data  was  continuously  acquhed  for  a  maximum  of  60 
seconds  after  the  initiation  of  dau  acquisition  to  capture  the  entire 
start  sequence.  The  total  pressure  and  temperature  measuremenu 
experienced  a  time  lag  of  tq>  to  one  second  behind  the  rotor  speed 
measurements,  but  this  did  not  affect  the  synchronization  of  mass 
flow  calculations  which  are  mdependent  of  rotor  speed. 


T-55  Rowpath  Geometry 
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ngur*  4.  T-55  DYNTECC  Flowpath  Model 
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Rgure3.  T55-L-712  Turboshaft  Engine 


THE  TC5-L-712  ENGINE 

The  T55-L-712  turboshaft  engine  (Fig.3)  is  a  gas  turbine  engine 
in  the  12  kg/t  class.  The  transonic  compressor  consists  of  seven  axial 
stages  and  one  centrifugal  stage.  It  uses  no  variable  geometry  but 
uses  a  single  start  bleed  over  the  sixth  stage  stator.  The  compressor 
operates  at  a  design  pennt  of  about  12  kg/i  mass  flow  with  a  pressure 
ndo  oi  about  8.0.  The  combustor  is  a  reverse  flow  annular 
combustor.  For  these  start  tests,  the  engine  power  turbine  was 
locked.  Hgure  4  shows  the  meridional  flowpath  of  the  engine  as  it 
was  implemented  in  the  DYNTECC  model 

Steady  state  instrumentation  included  total  pressure  and 
temperature  rakes  in  the  bellmouth  region  upstream  of  the  engine 
inlet  for  use  in  calculating  engine  mass  flow.  Pressure  and 
temperature  rakes  were  integrated  to  provide  overall  stage  pressure 
and  temperature  rise  characteristics.  High  response  instrumentation 
included,  but  was  not  limited  to,  a  string  of  high  response  pressure 


RESULTS 

Figure  5  shows  a  simulated  start  sequence  from  0%*80%  design 
corrected  speed  plotted  with  test  lesulu  from  a  successful  cold  start 
and  warm  start  Axis  scales  are  omitted  since  all  dau  has  been  non- 
dimensionalized  The  speedlines  in  this  plot  were  generated  from  the 
model  based  on  stage  characteristics  from  rig  test  data.  This  plot 
shows  overall  compressor  performance  from  compressor  inlet  to  the 
combustor  inlet  The  test  dau  goes  up  to  ground  idle  speed,  which  is 
approximately  61%  oi  design  corrected  speed.  The  two  experimental 
starts  follow  essentially  the  same  path  on  the  overall  plot  The  most 
obvious  feature  of  the  start  sequence  is  that  the  engine  appears  to  be 
operating  on  the  stall  line  through  most  of  the  start  The  simulated 
start  sequence  follows  a  slightly  lower  overall  operating  line  but 
closely  matches  the  ground  idle  point  The  effect  of  closing  the  start 
bleed  can  be  clearly  seen  in  the  simulation.  The  point  at  which  the 
exit  flow  becomes  choked  has  been  labeled  on  the  plot  to  emphasize 
the  smooth  transition  which  the  model  is  able  to  make.  This  is  a 
direct  result  of  the  pressure  /  area  boundary  condition  and  represents 
a  significant  enhancement  to  the  modeling  capabilities  of  DYNTECC. 
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Rgure  5.  Full  Start  Simulation  vs.  T-55  data. 
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Hgurs  6.  High  Rasponsa  Static  Prassura  Traca 


Figure  6  supports  the  observation  that  the  compressor  is  indeed 
operating  in  rotating  stall  during  much  of  the  start  sequence.  It  shows 
the  high  frequency  response  of  a  flush  mounted  shroud  pressure 
transducer  taken  when  the  engine  was  at  52%  of  design  corrected 
speed.  This  was  one  of  eight  circumf«entially  mounted  transducers 
at  the  same  axial  location  approximately  one  chord  length  upstream 
of  the  first  stage  rotor.  As  can  be  seen,  diis  sharp  pressure  rise  occurs 
with  a  fiequcncy  of  about  100  Hz,  roughly  60%  of  the  ground  idle 
compressor  spc<i  Further,  it  appears  apjwoximately  1.25  ms  later  in 
the  transducer  mounted  45®  in  ic  direction  of  rotation.  This  pattern 
continues  for  all  eight  of  the  transducers  in  the  array  in  front  of  the 
compressor.  It  seems  logical  to  assume  that  this  pressure  variation  is 
moving  ciroimfercntially  about  the  front  face.  Transducers  located  in 
an  axial  string  indicate  a  pressure  drop  behind  the  first  stage  that 
moves  with  this  pulse.  A  study  of  the  entire  transducer  pressure  trace 
indicates  this  pattern  exists  from  the  beginning  of  the  start  sequence 
with  up  to  four  of  tiiese  rotating  pressure  events  occurring  early  in  the 
start  sequence.  That  number  gradually  declines  until  the  event 
disappears  at  52%  ctf  design  speed.  This  apparently  rotating  pressure 
phenomenon  moves  circumferentially  at  between  40%  and  60%  of  the 
rotor  speed  throughout  the  sequence. 

Figures  7-12  show  the  start  sequence  on  a  stage-by-stage  basis 
for  stages  1-6  respectively.  The  background  specdlines  for  these 
plots  are  the  pre-stall  stage  characteristics  used  by  DYNTECX^  which 
were  derived  from  experimental  rig  test  data.  Axis  scales  have  been 
omitted  since  all  parameters  have  been  non-dimensionalized.  The 
experimental  data  for  stages  1-3  lies  on  or  near  the  stall  line  for  most 
of  the  start  sequence.  This  siqjports  the  observation  that  the  front 
stages  arc  operating  in  or  near  rotating  stall  throughout  most  of  the 
start  sequence.  The  model  closely  matches  this  behavior  and  matches 
the  ground  idle  point  quite  well  for  all  stages.  The  tendency  of  the 
model  to  follow  a  lower  operating  line  than  the  test  data  in  the  aft 
stages  is  not  clearly  understood  at  this  time.  Acrodjnamic  blockage 
is  inherent  in  the  stage  characteristics,  but  these  arc  only  valid  for  the 
bleed  schedule  which  was  used  during  data  acquisition  on  the  test  rig. 


It  has  been  postulated  that  an  explicit  treatment  flow  blockage, 
rather  than  accounting  for  blockage  in  the  stage  characteristics,  may 
improve  the  match. 

A  point  of  interest  is  how  the  ntodcl  is  able  to  run  at  speeds 
below  20%  of  design,  since  stage  characteristics  were  not  available 
below  20%  speed.  A  zero  speed  characteristic  was  input  to  the  model 
with  the  assumption  that  the  pressure  and  temperature  rise 
characteristics  would  be  zero  at  zero  speed.  The  model  uses  an 
interpolation  scheme  to  calculate  characteristics  at  intermediate 
speeds,  hence  it  interpolates  between  the  assumed  zero  speed 
condition  and  the  nearest  specdline  (20%)  for  the  very  low  speed 
calculations. 


CX3NCLUSIONS 

The  DYNamic  Turbine  Engine  Compressor  Code  (DYNTECC) 
siqjportcd  by  JDAPS  has  been  modified  to  incorporate  speed 
transients  with  pressure  /  area  boundary  conditions  allowing  dynamic 
simulation  ci  compressor  events  thiwghout  the  entire  operating 
regime.  The  enhanced  starting  model  has  been  configured  to 
simulate  starting  behavior  in  the  AUiedSignal  T55-L-712  axi- 
ccntrifugal  compressor  and  compares  favorably  with  start  test  data 
from  this  engine.  The  T55  compressor  normally  operates  in  rotating 
stall  until  it  reaches  approximately  50%  of  design  speed  (ground  idle 
is  about  61%).  This  is  near  the  starter  disengagement  point  for  the 
start  sequence.  During  the  start  sequence,  the  first  and  third  stages  of 
the  compressor  are  the  most  highly  loaded  axial  stages.  Suge  3  is 
most  highly  loaded  between  about  20%  and  43%  of  design  speed  and 
stage  1  is  loaded  more  highly  elsewhere.  Ownparison  of  interstage 
dau  with  the  model  shows  a  tendency  of  the  model  to  follow  a  lower 
operating  line  as  the  flow  approaches  the  start  bleed  location.  It  is 
believed  that  an  explicit  treatment  of  flow  blockage  may  improve 
these  results.  The  model  represents  a  significant  investigative  tool 
which  can  be  used  to  simulate  changes  in  start  bleed  schedules  and 
stage  matching  during  dynamic  events  “all  over  the  map”. 
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CORRECTED  FLOW 


FIgur*  7.  Stag*  1  Start  Saquenc* 


STAGE  2 


Figure  10.  Stag*  4  Start  Sequenc* 


Figures.  Stag* 2 Start S«qu*nc* 


Figure  9.  Stage  3  Start  Sequenc* 


Figure  11.  Stag*  5  Start  Sequenc* 
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